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Abstract 

We consider the small velocity sum rules for heavy flavour semileptonic 
transitions that are used to estimate the zero recoil values of semileptonic 
heavy flavour form factors. We analyze the complete 0 ( 015 ) radiative cor¬ 
rection to these sum rules. The corrections are universal and influence all 
’’model-independent” bounds previously derived for semileptonic form factors 
at zero recoil. 
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1 .Introduction. 

An accurate determination of the Kobayashi-Maskawa matrix element V cb is one of 
the most important tasks of the heavy quark theory in the physics of heavy hadrons. 
As is well known there are at least two possibilities for measuring i) to measure 
inclusive semileptonic decays of heavy hadrons or ii) to extrapolate differential decay 
distributions of exclusive semileptonic transition to the zero recoil point that gives 
us | V c b \ Ja- In the second method precise theoretical predictions for the values 
of form factor / 4 at zero recoil are extremely important for extracting V c b from 
experiments. 

A general approach for obtaining such predictions was suggested by Shifman et al. 
[1,2]. It is based on recent progress in the analysis of inclusive semileptonic decays of 
hadrons containing one heavy quark, where the operator product expansion (OPE) 
method and HQET were applied [3-7]. The leading order result agrees with the free 
b-quark decay picture. Nonperturbative corrections start to appear only at order 
( A a/q° ) 2 an d are determined by the matrix elements of only a few local operators 
such as the operators of the chromo-magnetic and the kinetic energy [3-7]. 

The procedure of Shifman et al. [1,2] consists in considering moments of spectral 
distributions in the small velocity (SV) limit which allows one to obtain corrections 
to known sum rules as well as new sum rules. A very important result of this 
approach is that now one has an estimate of the deviation of the value of the exclusive 
form factor / 4 from unity at zero recoil. 

The aim of the present note is to analyze the 0(as) corrections to the SV sum 
rules. Note that part of these corrections had already been incorporated in the 
original derivation of the sum rules. We are referring to the vertex renormalization 
factors r/Ay which correspond to the finite renormalization of the vertices (c|c 7 Al 5 | 5 ) 
and (c|c 7 m 7 5 6 | 6 ), respectively. To the best of our knowledge these factors were first 
introduced in [3]. However, we would like to emphasize that the vertex renormal¬ 
ization factors are not the whole story at the 0(as ) level and extra care is needed 
to derive correct sum rules with 0(«s) accuracy. The essential point here is that 
the diagrams with two partons in the intermediate state are also involved. 

2. Derivation of SV sum rules. 

Let us begin with the forward scattering amplitude: 

T^iqv) = -i J dxe~ iqx < H h \ j^x)*j u ( 0) | H b > (1) 

Here j v is the appropriate current and v is the velocity of the heavy hadron H b . 
The function T(qv ) is an analytic function in the (qv)— complex plane with the 
appropriate cuts. The structure of the T IIU cuts is shown in Fig.l. The cut 0 < 
qv < Mu b — Mh c corresponds to the decay of heavy hadron H b —> H c + ul while 
the lower cut qv < 0 and the upper cut qv > Mu b + Mh c represent other (crossed) 
physical processes. Following Ref.[l] we shall argue latter on that, using duality 
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concepts, the contributions from the latter ” crossed” cuts to the sum rules can be 
neglected in as much as the theoretical and phenomenological contributions on these 
cuts can be equated to one another. The imaginary part of the forward scattering 
amplitude on the ”physical cut” 0 < qv < M Hh — M Hc determines the hadronic 
tensor W IJtv = — 4 ImT IJV which in turn determines the inclusive decay width of the 
hadron. 

Before discussing the calculation of the 0(a s ) corrections we would like to men¬ 
tion that there exist two different approaches for deriving SV sum rules. One was 
proposed by Chay, Georgi and Grinstein (CGG approach) [4] and the other by Shif- 
man et. al. (BSVU approach) [1], Both approaches are based on the duality idea 
(global and local) but use a somewhat different language. We shall discuss them in 
turn. 

The basis of the CGG approach is connected with the possibility to perform 
an analytic continuation of the forward scattering amplitude to the whole complex 
(qv )—plane and to connect the integral over the physical cut with the integral over 
a ” large-radius” circle in the complex plane where the OPE is justified. This can be 
thought of as a formal statement of the assumed duality. A representative integration 
path Ci is shown in Fig.l . 

Integrating over the physical cut and equating moments of theoretical spectral 
functions with their phenomenological counterparts we obtain the sum rule: 

[ de(e) n W phen (e) = f de(e) n W QCD (e), (2) 

J phys.cut J 

where e — qv — M Hb — M # c . The integrand on the r.h.s of these sum rules includes 
nonperturbative ^ power corrections as well as radiative QCD corrections which 
can be systematically incorporated by using standard OPE calculations. 

In the second (BSVU) approach [1] one writes down a dispersion relation repre¬ 
sentation for the forward scattering amplitude 


TM = i / (3) 

and expands the integrand in powers of where, in the rest frame of the initial 
hadron, e = mb — m c — q 0 . Such an expansion is justified only for e-values e » Aq C d 
and e « 2 m c . However, in general there are excited states with large invariant 
mass (~ m c , (nib — rn c )) whose contributions are not suppressed (the perturbative 
analogue is surely ’’hard” gluon emission). This means that one cannot expand the 
integrand of Eq.(3) in terms of powers of It is then inevitable to split the region 
of integration in the dispersion integral at some scale /a >> A qcd and /a << 2 m c . 
Rewriting Eq.(3) as 





r d m 

J o e — e 


fM Hb -M Hc W(e) 
+ / de — ±4. 

J p 6 — 6 


(4) 
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and assuming e >> /x, we can expand the first term on the r.h.s. in terms of powers 
of / while the second term on the r.h.s. (which reflects the contribution of excited 
states) does not generally have a Laurant expansion for e << /i. However, at this 
point one can invoke local duality to evaluate the second term in the r.h.s.. In some 
sense this piece of the theoretical spectral function has to incorporate all radiative 
corrections coming from the hard region. In other words we apply duality to equate 
the integral over the region [/x, Mfj h — Mh c ] on the partonic side to the sum over 
excited states (with energy E > /x). Thus we assume that duality is valid up to the 
scale n « m c . Then doing an expansion in e we obtain the result: 

^ de(e) n W Phen -(e) = dl(e) n W QCD (e). (5) 

J o Jo 

We want to emphasize that /x is the scale where the excited states come into play. 
In this sense the/x scale is quite similar to the energy of the continuum threshold 
in the usual QCD/SVZ sum rules. It is well known that a particular choice for the 
onset of the continuum contributions affects the final result of the QCD/SVZ sum 
rules analysis and typically results in 20-30 percent uncertainties. This shows that 
some care is also needed to estimate the size of this uncertainty in the analysis of 
the SV sum rules. 

Remember that, when calculating QCD radiative correction to the Wilson coeffi¬ 
cients, we have to introduce yet another scale which divides the region of integration 
into a perturbative and a nonperturbative region. In order to distinguish these scales 
we shall denote latter scale as fioPE whereas we denote the duality scale discussed 
above by fi D . From perturbative QCD we roughly expect /Jope > 300MeV whereas 
we take /X£> ~ 1 GeV for the duality scale. By definition the sum rules (2)-(5) cannot 
depend on the scale i^ope- Also the SV sum rules do not depend on /X£> to leading 
order in as and to any arbitrary order in — because the spectral functions originate 
entirely from the resonance region. The leading order spectral function (leading in 
as) consists of a 5(e ) function and 5-function derivatives. Note though that, in the 
next to leading order in as, the spectral function starts to depend on fin explicitly. 

Thus it is impossible to derive realistic ’’model-independent” bounds for the form 
factors: the choice of ’’switching on” the /xj>— dependence is intrinsically model 
dependent. It is worth mentioning that the same problem appears in the CGG 
approach and is, technically, connected with the possibility to choose different inte¬ 
gration contours in the (gv)-plane, each contour C -2 being defined by the point where 
it leaves the physical cut - this is the aforementioned jiu ambiguity (see Fig.l) in a 
different guise. 

It is important to realize that the dependence on /jd appears only at next to 
leading order in as, and not at leading order as in the usual QCD/SVZ sum rules. 
To estimate the size of this ambiguity and thereby to understand the accuracy of 
the SV sum rules we have to compute QCD radiative corrections to the SV sum 
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rules at arbitrary values of the scale /i/> 

3. Results for structure functions. 

In this section we present the result of calculating QCD radiative corrections to those 
spectral functions that are needed for the zero velocity sum rules. The hadronic 
tensor W^ v can be expanded in terms of 14 structure functions (see [4,6,8]): 

W^ = -g^Wi + v fl v u W 2 - ie llua pv a q !3 W 3 + q^W* + (q»v v + q u v fl )W 5 (6) 
-qs[~9^ W 6 + Vfj,v u W 7 - ie^ a pV a q l3 W% + q^Wg + (qypy + q^v^Wio] 

+ SvvjWu + {s^q u + s u q^ L )\\ i -2 + i^^uapWs^Wi^ + ie fllia pq a s^Wi4 ! 


where v is the velocity and s is the spin of the initial hadron. 

As a next step one defines diagonal hclicity structure functions Wl (longitudinal), 
Wt lr (transverse left,right), and Wo (time-component or scalar) by contracting 
the hadronic tensor with the appropriate polarization vectors n u x n u * x (A = 
0, +, (see Ref.[ 8 ] for details ). 

At zero recoil one finds: 

W L = Wi 

W 0 = -W 1 + W 2 + (qv) 2 W 4 + 2(qv)W 5 (7) 

Wt L zR = Wi ± nsW TS , (8) 


where 

Wts = (W 13 + qvWu) (9) 

and n defines the quantization axis of the off-shell IT-boson. 

Next we calculate the zero recoil O(as) contributions to the hclicity structure 
functions ITo, Wl and Wts- The generic graphs which contribute to the 0(a s ) 
correction are shown in Fig.2. There are in principal several possibilities to cut 
the graphs in Fig.2. The difference is the number of partons in the intermediate 
state. The cuts with one parton in the intermediate state reproduce the vertex 
renormalization t]a,v factors which were discussed earlier. The two-parton cuts give 
rise to the absorptive contributions in the range 0 < qv < mj — m c . We mention 
that the infrared singularites cancel in the sum of the two-parton intermediate state 
contributions in the soft gluon limit qv —> rrii, — m c . This is in accord with the 
observation that the virtual corrections have no infrared singularity at the zero 
recoil point [3,10]. 

We first consider the correlator of two axial vector currents = 07^756 and 
obtain the following zero recoil contributions: 


W L AA (t) 


v 2 A wE orn (t) + 


asCp (5 1 2 + 10 tx + 3x 2 + 2x + 3) (t + 2 x)t 
2n 6 (t + x) 3 rrib 


( 10 ) 
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( 11 ) 


W 0 AA (t ) = rj 2 A W 0 Born (t) + 
W TS AA (t ) = ^W^i) + ■ 


(XsCf (3 1 2 + 6 tx + x 2 + 2x + 1 )(t + 2x)t 
2n 2(t + x) 3 rrib 

&sCf ( t 2 + 2te + x 2 — 2x — 3)(t + 2x)t 
2n 6 (t + x) 3 rrib 


( 12 ) 


where t = — m , c — and x — — . The vertex renormalization factor in Eq.(12) reads 

mb mb 1 v / 

[3]: 


VA = 1 - — (j^log(x) + (13) 

7T 1 — X 3 

and, for the corresponding vector current case, 

W = 1 - + 2). (14) 

7T 1 — X 


The 0 ( 0 ^ 5 ) Born contributions (which are state-dependent) include both the zeroth 
order result and the effects of non-perturbative corrections and can be found else¬ 
where [1,2,6-9]. Concerning the vector current case the structure functions (10)-(12) 
are the same except for a sign change in the fourth term of the Erst brackets (2a;- 
term) in all formulas (10-12). 


4. QCD Radiative corrections to SV sum rules. 

Next let us discuss the 0(a s ) corrections to the SV sum rules. As concerns sum 
rule applications the three helicity structure functions Wl AA , Wts AA and Wq 1 
are the most important. These are the structure functions that have nonvanishing 
contributions at zero recoil from the quasi-elastic contributions B —> D, D* and 
A b —> A c and f \ ^ We want to discuss how these sum rules are modified 

when 0(a s ) radiative corrections are taken into account. 

For the ” time-component” helicity structure function IT 7 ,} 1 we obtain 

(/ +/ )Wq A (t)dt = i] V + — —(Ji(/xd)- Ji(noPE)) + n/pert(no PE ), (15) 

J 0 ''HOPE 

where Ji is defined in the Appendix. The last term in Eq.(15) is a symbolic no¬ 
tation for the non-perturbative parameterization of the contribution from the soft 
region: [0 , Hope]- Since we are mainly interested in how previously derived sum 
rules change when the new 0(a s ) corrections are included the last term in Eq.(15) 
will not be written out explicitly in the following since these contributions have been 
investigated in previous papers (see Refs[l,2,8,9]). On the other hand, Jq d Wq VV dx 
is connected with the sum rule for the vector form factors, e.g. in the Af, case (for 
form factor definitions see e.g. [8]): 

I E/M 2 < rwv^dx (16) 

4 = 1 J ° 
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Let us just for illustrative purposes assume hope —> 0 and Ho = m>b — m c that 
corresponds to integration over whole physical cut. We then obtain: 

|£/v| 2 < Vv + as ^ F • ( log(x){hx 2 + 2x - 1) - (x 4 + x 3 - x - 1)). (17) 

o7T 

But integrating over resonance region hd << m c and assuming hope —>• 0 we have 

I /ul 2 < Vv + - l) 2 -^f + n/pert(0). (18) 

^ 4tt m2 

Now let us consider the SV sum rules more carefully and estimate their dependence 
on the duality scale hd and Hope■ The basic function Ji(h) is shown in the Fig.3. 
Let us emphasize that J\ is small in the infrared region and hence our result for 
perturbative correction is not sensitive to Hope- Thus we may set Hope = 0 and 
forget about the Ho p E-dependence of operators with higher dimensions and use the 
known expressions for the r)v,A factors [3]. 

In order to make reliable estimates we have to decide on the value of hd- hi 
principle using the CGG approach we may choose any hd in the interval [hope, mb — 
m c ] based on various assumptions about the applicability of duality. Conventionally 
one takes hd — 1 — 3 GeV in QCD sum rule applications. We see from Fig.3 that the 
result for the QCD radiative corrections to leading operator does have a substantial 
dependence on hd in the region 1-3 GeV. The size of the effect varies from 0.5 % 
at hd = lGeV to 3. % at hd — 3GeV, here we have used a s = 0.3, mb = 4.8Gev , 
m c = 1.5GeV for definiteness. 

A similar situation occurs for all the other sum rules considered in the literature 
before. For instance, taking the external current to be axial and projecting on the 
longitudinal helicity function for B D* , the radiative corrections considered here 
change the original inequality | f A | 2 <| r] A | 2 +n/pert derived in Refs. [1,2] to 

|/a| 2 < Va H—^— ■ (J 2 (pd) — J 2 {pope )) + n/pert{noPE) (19) 

Ztt 

with J 2 from the Appendix. In the case Hope 0 and hd « m c we obtain 

\Ja\ 2 < Va + —-(^ 2 + \x + 1)^§ + n/pert{ 0). (20) 

47T 3 m z c 

This result changes the prediction for the bound f A < 0.94 obtained in Ref. [2], The 
bound is pushed upward by 0.5% and 2.% for Hd — 1 and 3 GeV, respectively. 

Note that the authors of Refs. [1,2] have used one single scale h = Vd = Pope 
implicitly including considered effects to nonperturbative operators. But if this is 
the case one has to estimate and include the /^-dependence of the nonperturbative 
matrix elements. This is important at 0(a s ) accuracy since the //-dependent part of 
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these operators turns out to be of the order as( ^-) 2 ~ cks at // m m c . A simple way 
to avoid this problem is to take /i « m c as it has been done in Ref. [1,2]. However, 
strictly speaking, this choice for the value of /j is not quite harmless as one applies 
duality concept up to an extremely small scale // << m c where local duality can 
break down (see for example Ref. [11]). Moreover, if we use the commonly accepted 
numerical values for A qcd and m c then the reliability of the strong inequality is 
doubtful: 

A qcd << h << ( 21 ) 

On the other hand a numerical analysis of our formulas shows that the typical 
size of radiative corrections to the small velocity sum rules is of the order of 1 % at 
H = m c . We can regard this as an estimate of the theoretical uncertainty within the 
SV sum rule calculations. Phenomenologically this results in a 1 — 2% increase of 
the relevant bounds on the zero recoil form factors. 

Applying SV sum rules techniques one can also obtain inequalities for the matrix 
elements of non-perturbative operators. For example for the A^-case discussed in 
Ref. [ 8 ] the following inequality was obtained: 

R, 2 + ^ < 0 (22) 

Here /i n 2 stands for the expectation value of the heavy quark kinetic energy oper¬ 
ator while n s 2 parameterizes the Aj correction for the axial vector current matrix 
element between Af, states. 

Using our previous results the inequality (21) gets changed to 


Rs 2 + - ^ F '^(^(Ad) — J^ope)) < 0 

O 27T 


(23) 


when the 0(as) radiative corrections are included. The function J 3 is defined in the 
Appendix. Again at hope —► 0 and hd « m c we obtain 


^ ^ + ^(1 + l x - i A < 0 . 


(24) 


All terms in the last equation can be of the same order of magnitude in principle 
(fj, 2 ~ 0.5 GeV, last term in Eq.(24) is about 0.8(^) 2 ) and there are no a priori 
reason to discriminate between them. But then the values for /i 2 and fi n 2 are not 
connected directly what is the major advantage when the radiative corrections are 
neglected. 

5. Conclusion 

To conclude, we have analysed the 0(a s ) radiative corrections to zero recoil sum 
rules for semileptonic heavy hadrons form factors. Our results are universal and 
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shift all previously derived model-independent bounds on zero recoil form factors by 
1-2 percents with obvious consequences for the extraction of V c b from experiment. 
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clarifying conversations on the subject. 

7. Appendix 

In this appendix we collect our results for the 0(as ) QCD corrections to the lead¬ 
ing order term in the SV sum rules. These corrections arise from the two parton 
intermediate states. 

The QCD radiative correction to the zeroth moment of the Wo 1 1 hclicity structure 
function ( J\ (/i) = Jq Wq V ( t)dt) in the vector current correlator is 

JM = T 7 —~(— 2 (u x + log(^^-){y + a;) 2 )(5a: 2 + 2a: - 1) (25) 

4 [y + x) z x 

-y 2 (l7x 2 + 2x - 1) - 3 y 4 - 12 y 3 x) 

with y = y/rrib and x = For the W^ AA helicity structure function the correction 
to zeroth moment is given by: 

Mv) = TR7 ~r w (~ 2 (y x + + x) 2 ){7x 2 - 2x - 3) (26) 

1 2{y + x) z x 

-y 2 (27x 2 - 2x - 3) - 5 y A - 20 y 3 x). 

Finally for the transverse structure function Wts AA one has 

JM = 19 / \ s 2 (~ 2 (y x + l °9(^^){y + x) 2 ){x 2 + 2a: + 3) (27) 

12 (y + x) z x 

—y 2 ( 5x 2 + 2a: + 3) — y 4 — 4 y 3 x). 
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Figure captions 


Fig.l The structure of the forward scattering cuts in the qv complex plane. C\ 
and C '2 are representative integration paths. 

Fig.2 The generic graphs which contribute to the 0(as ) correction. 

Fig.3 Dependence of Ji(fj) (upper), J 2 (/x) and (lower on r.h.s.) on /j at 

rrii, = 4.8 GeV, m c = 1.5 GeV. 
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